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Abstract 

So far, it is not well known how to deal with dissipative systems. There are many paths of 
investigation in the literature and none of them present a systematic and general procedure to 
tackle the problem. On the other hand, it is well known that the fractional formalism is a powerful 
alternative when treating dissipative problems. In this paper we propose a detailed way of attacking 
the issue using the fractional calculus to construct an extension for the Dirac brackets in order 
to furnish the quantization of nonconservative theories through the standard canonical way. We 
believe that it can be the first step to construct gauge theories from second-class nonlinear systems 
using these extended Dirac brackets. 
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I. INTRODUCTION 

Very popular in the nineties, where an industrial production of papers concerning methods 
treating constrained systems, the Dirac brackets (DB) pj] were an unmodified common point 
between all papers in the subject. The main objective of many works were to convert second- 
class systems in a first-class one, which is considered a gauge theory, i. e., the holy grail for 
the Standard Model. Although not so popular as before, the analysis of constrained systems 
deserves some recent attentions in the literature j2[. 

In few words we can say [1] that the main feature of gauge theories is the existence 
of constraints which fix boundaries in the phase space of gauge invariant systems to a 
submanifold. In fi\ Dirac covered all the main issues concerning constraint systems. Namely, 
a Hamiltonian approach to gauge theories and general constrained and, consequently, the 
corresponding operator quantization frameworks. Later on, the path integral method was 
found to be useful for quantizing gauge theories j3] and the so-called second-class systems 
[HI, where the conventional Poisson bracket must be replaced by the Dirac bracket in the 
quatization procedure. 

However, in constrained systems, it is possible to solve constraint equations [3|. The 
formalism proposed by Dirac for classical second-class constrained systems uses the Dirac 
bracket to deal with the evolution problem. The procedure is to apply the Dirac bracket 
to functions of canonical variables in the unconstrained phase space, which avoid problems 
concerning the restriction of systems to constrints submanifolds [sj. 

On the other hand, there are various problems when considering classical systems besides 
the ones involving the quantization of second-class systems as we saw just above. These 
problems constitutes the so-called nonconservative systems. The curiosity about them is 
that the great majority of classical systems is nonconservative and nevertheless, the most 
advanced formalisms of classical mechanics deals only with conservative systems 



Dissipation for example, is present even at the microscopic level. There is dissipation in 
every non-equilibrium or fiuctuating process, including dissipative tunneling 0], electromag- 
netic cavity radiation jsf and so on p]. 

Through the years, the number of methods concerning nonconservative systems. However, 
none of them show the same effectiveness found in the mechanics of conservative systems. 

One way to attack nonconservative systems is through the FC since it can be shown that, 
for example, a friction force has its form resulting from a Lagrangian containing a term 
proportional to the fractional derivative which is a derivative of any non- integer order . 

Fractional calculus is one of the generalizations of the classical calculus. It has been 
used in several fields of science such as turbulence and fiuid dynamics, stochastic dynamical 
systems and plasma physics |9i]. In this work we will use the well known FC to analyze the 
well established DB. The objective is to construct a generalized DB capable of treating a 
bigger number of mechanical systems than the standard DB. 

Since we believe that the FC is not well known yet, we tried to construct a self-sustained 
paper so-that the issues are distributed as follows. In section II we furnish a short his- 
tory about FC together with its main equations and formulations. We will follow here the 
Riemann-Liouville (RL) approach. In section HI we establish the so-called fractional varia- 
tional principle, the ground stone for our cherished result. However, we have to perform a 
modification of this fractional principle in order to include constrained systems. This was 
accomplished in section IV. Next, we establish different initial conditions to obtain a general 
formulation for the Dirac description for constrained systems. In section V, we use the free 
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relativistic model to apply the fractional bracket. As ususal, the conclusions, perspectives 
and last comments are depicted in the last section, the sixth one. 



II. FRACTIONAL VARIATIONAL PRINCIPLE 

It is well known from the current literature that the fractional approach can describe 
more precisely a myriad of physical systems and that the formalism can be incorporated 
easily in many classical and quantum systems. Also, its use can be extended up to field 
theory domain. Although the fractional calculus technique is still not very popular, it can 
be useful to study several problems from different areas (besides the ones mentioned last 
section) of science such as viscoelasticity and damping, glassy condensation, diffusion and 
wave propagation, electromagnetism, chaos and fractals, heat transfer, biology, electronics, 
signal processing, robotics, system identification, genetic algorithms, percolation, modeling 
and identification, telecommunications, chemistry, irreversibility, control systems as well as 
engineering, economy and finance 11 23l . 



The generalization of the concept of derivative with non-integer values goes back to the 
beginning of the theory of differential calculus. Nevertheless, the development of the theory 
of FC is due to contributions of many mathematicians such as Euler, Liouville, Riemann, 



and Letnikov M26 



Since 1931 when Bauer [27| showed that we can not use the variational principle to obtain 
a single linear dissipative equation of motion with constant coefficients, a new horizon of 
possibilities were glimpsed. Nowadays it has been observed that in physics and mathematics 
the methodology necessary to understand new questions has changed towards more compact 
notations and powerful nonlinear and qualitative methods. Derivatives and integrals of 
fractional order have been used to understand many physical applications. For instance, 
questions about viscoelasticity and diffusion process may have a more detailed description 
when this approach is used. In nature, the majority of systems contains an internal damping 
process and the traditional approach based on energy aspects can not be used everywhere 
to obtain the right equations of motion. 



So, after Bauer's corollary, Bateman [28[ proposed a procedure where multiple equations 



were obtained through a Lagrangian. Riewe [29[ observed that using FC it was possible to 



obtain a formalism which could be used to describe both conservative and nonconservative 
systems. Namely, using this approach one can obtain the Lagrangian and Hamiltonian 
equations of motion also for nonconservative systems. In jl^] Agrawal studied a fractional 
variatonal problem. A fractal concept applied to quantum physics has been investigated and 
reported in (30j . 

Recently this subject was revisited in jlH and the solution of a fractional Dirac equation 



(order 2/3) was introduced in [32 



Modified Equations: To begin with, let us consider the action functional below, defined 



by means of RL fractional derivative in a configuration space [33 



S[X] = w^ / mr),q{r),r){t-rr-'dT, (1) 

where T{ai) is the traditional Euler gamma function, with a G (0, 1) and g = ^ is the 
derivative with respect to the intrinsic time r G {a,t') and t e [to,t'] is the time for some 
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observer in a particular referential. Obviously when a — )■ 1 we re obtain the usual functional 

Six] = jmr),q{r),r)dT. (2) 

If the curve x is an extremal of S[x\ the sufficient and necessary condition required is 

dL d f dL\ f I — a\ dL „ . , 

— = 0, 2 = l---n . (3) 



dqi dr \ dqi J \t — t J dqi 

The Euler-Lagrange equation above, for some action fractional functional must be obeyed. 
However, if we consider now the same action variation in phase space 

6S = 0, (4) 

we have that the integral 



6S = - 



^5 [[pq- H{p„ g„ t)] (t - rr-'dr = (5) 



will permit us to write a new set of perturbed equations. 



1^ = IT- (6) 

dpi 

dH fl-a\ _ 

Pi = + Pi z ' \V 



dqi \t-T 

which can be understood as the Hamilton- Jacobi equations when this new action functional 
is considered. It is clear that when a — )■ 1 our results will return to the usual shown 
above. 

The presence of a fractional factor 

1 — a 

is responsible for the generation of a time-dependent damping into the dynamics of the 
system, which is very useful to study models with smooth turbulence. Furthermore it is 
possible to establish a relationship between the fractional Rayleigh dissipation function {R) 
and the Euler-Lagrange equation [33] . 

dL _ d fdL\ _ dR _ ^ ^_ _^ ^ 
dqi dr \dqi) dqi 



where 

R = l( 

t - T 



R = L[^-^]. (9) 



Note that in ([8]), the dissipation function is part of the extended Euler-Lagrange equation. 
However, the origin of the third term is non-standard due to fractional analysis. 
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III. MODIFIED VARIATIONAL PRINCIPLE ON CONSTRAINED SYSTEMS 



Now our main objective is to obtain an extended analysis which allow the quantization 
of classical systems with turbulence flow in field theory. We know that the quantization of 
a classical field theory in a natural context is not a straightforward unique process. The 
replacement of classical Poisson brackets by commutators of quantum operators can not be 
carried out simultaneously for all conceivable dynamical variables without paying the price, 
i. e., internal obstructions will occur 36, 37l . 



In general, the commutation formalism is restricted firstly to a certain class of variables, 
such as the canonical coordinates of the theory. All commutators obtained will be derived 
from this first set. However the classical theory may be substantiated in terms of any set of 
canonically conjugated variables in such manner that the transition from Poisson brackets 
to quantum commutators leads to a "weird" quantum theory, depending on the chosen 
canonical coordinates system. 

This kind of problem usually occurs when the classical theory has constraints, and the 
right prescription for this was first formulated by Dirac and Bergmann 36|, l37|], where they 
pinpointed the right bracket algebra to be used. Thus our goal now is to extend our last 
result to constrained systems. The action then can be considered in phase space. 



S 



1 



pq 



(10) 



where H is 



H = H + X„ 



(11) 



The question involved in such systems is that when we carried out the Legendre transfor- 
mation (mapping) where, 

L{q,q,t)^H{p,q,t), (12) 

and to define the canonical momenta as = |4 perhaps the N quantities are not all 
independent functions of the velocities. We can not eliminate the g^'s and obtain 



(pa{q,p) = , a 



M 



(13) 



which are the constraints equations. 

Extending our discussion then, we write the variation for ( JTOl) as 



55" 
1 



1 

dH 



pq - H{pi, qi, (f)a) {t - r)" dr 



SpiQ - {Pi-Pi{ -j—^^ ^^i^ - 



6pi + Xa^^^Qi - Xa^T—5pi 

opi dqi dpi 



{t - r)"-Mr 







(14) 
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After some algebraic manipulations some terms can be isolated allowing us to write the 
Hamilton- Jacobi equations for the fractional constrained case, 

. _dH d<Pa 

Qi — -7; r Afl-;- — , 

opi dpi 



dH 



Pi 



-Xa- 



+ Pi 



a 



t - T 



(15) 



% dqi 

and again we have in second equation of f|T5|) a second-term representing the fractional 
contribution. 

Dirac Bracket: Consider a dynamical variable Q[qi,Pi,t] and using (|T5l) is obvious that 



de dQ de dQ 



dqk V dp^ dp^ 



+ 



dQ 

dpk 



dH 

dq^ 
dQ d 



X 



dq^ 
dQ 



a. , k 



a 



t - r 



ae _ d^dH_\ 

dqk dp'' dpk dp'' ) ^\ dqt dp'' dpk dp'' 



Pk 



{a - 1) dQ 
[t — r) dp'' 



+ 



dQ_ 
'dt 



{Q,H] + \a{QAa]-Pk 



1) dQ dQ 
'dt 



{t — r) dp'' 

The constraints are dynamical variables too. Then, substituting some of the constraints in 
( !T6|) we have that. 



dt 



{0a, H} + Xb{4>a, (f)b} - Pk 



[a-1) 



+ 



{t — r) dp'' dt 



(17) 



and solving for Aa, we finally obtain a new result for the DB, on a fractional context, namely, 
{F,Gr = - {F,<Pa}C;,' G} + 



{a - 1) d(f)t 
{t - r) dpk 



Pk- 



{a - 1) dF 
{t - r) dpk ' 



(18) 



Our calculations show precisely this new result as a natural extension for the DB. We must 
observe that the usual DB appears inside the fractional correction, the matrix Cab = {4'a, 4'b} 
is the constraint matrix. It is obvious that, when a — )■ 1 we re-obtain the usual approach. 



IV. FRACTIONAL EMBEDDING 



Our next step is to build a general way to obtain the Dirac description for constrained 
systems. For this we will consider the problem under different initial conditions. A different 
and more general approach to analyze any dynamical system begins by considering the 
action as a function of generalized coordinates [29| . 

S[x] = J^L{q:{T),Ql,{T),T)dT (19) 

q: = {aDtrXrit) , QI,{t) = {tD^r'xrit), 
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with r = 1,2, ... R coordinates considered, n = 1. 2. ... M is the sequential order of the 
derivatives for the generahzed coordinates q and n' = 1,2,... M' is the same for the coor- 
dinates Q. It can be showed that the necessary condition for an extremum of S is satisfied 

by 

i^|:(.-n-f .W^^o. (.0) 

and the momenta have the following form 



Pn- 2^ \t^h ) r. > 

k=n+l 

-TT-r ( r\o.\k-n'-\ 



(21) 



k=n'+l 



It is important to observe that we could extend the approach to a phase space just 
considering the usual action functional depending on the generalized fractional coordinates. 

The Dirac formalism can be easily obtained here. It is well known that it is useful in 
Lagrangian constrained systems. Now we propose its extension using the fractional calculus 
to encompass constrained non-conservative systems. Of course we could define our initial 
conditions in a different way and consequently to obtain other final expressions. We real- 
ize that it is a very general form to deal probably non-linear systems and other kinds of 
phenomena. With this objective we define our constrained Hamiltonian, 

^ = i^+^A,efe + ^«, (22) 

k k' 

where now the constraints are in fractional form also. Namely, they can be expressed in 
terms of fractional formalism. We will define them by means of the RL prescription, 

(23) 

The resulting action is, 

/t' R N-l R N'-l 

r=l n=0 r=l n'=0 

(24) 




and using of the variational principle again. 



5S = 0, 



(25) 
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we can calculate the Hamilton- Jacobi equations, 



dH , 9$ 



k 



OH dXy 

These forms for the Hamilton- Jacobi equations are new in the literature and introduces an 
extension of the Poisson bracket into the RL context presented in [38i] . It is natural that the 
next step is to obtain the proper DB expression. One way to do that is to consider some 
dynamical variable F{q'^,p^,Q^,,nl') where 



dF BF BF BF BF BF 

and after using ( 126|) it is straightforward to build our final and main result for the DB in 
RL context, namely, 

{A, BY = {A,B} - {A,cP,}C^,'{<j)i,B} - {A,Xk'}E^,l{xi',B} (28) 
where C and E are constraint matrices as in the standard Dirac constraint formalism. The 
consequently quantization can be described also in the standard way as 

[A,B] = ih{A,By . (29) 

From the moment we have constructed a proper fractional form for the DB, we believe 
that the conversion methods for obtaining first-class systems from second-class ones for 
non-linear models is a consequence. 



V. THE RELATIVISTIC FREE PARTICLE 

Our objective in this section is to show an application (new in the literature) of the 
fractional calculus and after that, of the fractional Dirac bracket introduced in fl28|) . 

To fix our ideas let us consider a simple example and we will use the relativistic free 
particle model to apply the fractional embedding. This model is well known, and its usual 
Lagrangian is given by, 

L = — (30) 
Using the ideas of the last sections, the action under consideration is, 

L = -m^{{aD-ty-Xrf + {{tD-,y-'xrf. (31) 

Geometrically speaking, fractional models provide us with a memory effect in convolution 
integrals and give us some differential equations with bigger expressive power. This allow 
us the consideration of several different physical situations as viscoelasticity as well as more 
abstract scenarios as mapping using tensorial fields. Physically we can understand the 
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{1 — a) — th Riemman- Liouville derivative order of some individual velocity v as the same 
velocity Vob from the view point of the independent observer 

We are going to restrict our calculations for the case when a = = n' = 1 therefore 
we have now, 

L = -m^ {{aDK)xry + i{tDh)xry (32) 

or in a simple way, 

L = -msJ{q[Y + {Q\Y- (33) 

Of course we could consider different orders to the derivative operator, but our main inten- 
tion now is to apply the method and show its usefulness to quantize canonically fractional 
systems. 

Using the definition of generalized momentum we obtain the two conjugated momenta. 
The primary constraint is, 

01 = pl + nl-m^ ^0 . (35) 

As the canonical Hamiltonian is zero, to construct the extended Hamiltonian by Dirac's 
prescription we can write, 

H = \(pl + nl-m^), (36) 
and since it is a first class constraint 

01 = {(PuH} = 0, 

there is a trace of gauge symmetry, therefore we do need to fix the gauge and our choice is 

g? + Q?-r = 0. (37) 
Now we have two second class constraints 

{01,02} = -2(p° + 7r°), (38) 
and the new extended Hamiltonian can be read as 

H = + {<? - + HWif + {Q\f - r) (39) 

the time evolution of these constraints give us the right form for the Lagrange multipliers, 
and the extended Hamiltonian in its final form is. 

Using the definition of Dirac brackets we can write finally that, 

{Q\,nir = 6^^ - (41) 



Fractional Dirac Bracket and Quantization for Constrained Systems 



10 



and the quantization is directly obtained through (129|) . 

Finally, we have to say that it is a first paper in the literature using fractional calculus 
to investigate nonconservative physical systems. It is obvious, in this first paper condition, 
that we are beginning to understand the physics behind the fractional formalism. To help 
us in this task, the next step would be to study a solid nonconservative system like the 
radiation damping, which is not understood using the standard (non-fractional) formalism. 
It is the target of current research by the authors and will be published elsewhere. 

VI. CONCLUSIONS AND PERSPECTIVES 

In conclusion, we proposed two kinds of fractional formulation for Dirac brackets. The 
first one was based on RL definition but incorporated directly to the action functional. We 
showed that Hamilton-Jacobi equations appeared deformed by the fractional contribution. 
Consequently the Dirac bracket also suffers the same kind of modification. However this 
first approach does not seems to be the right way to treat such questions. Subsequently 
we changed the formalism considering then an usual form to the action but redefining the 
coordinates now in a generalized prescription using the fractional definition as Riewe's pre- 
scription. The constraints were defined in the same way and the consequence was the 
extension of the usual Euler-Lagrange equations of motion to the fractional scenario. 

Finally, we obtained the final form for the fractional DB which has an additional term 
due to the FC contribution. The standard DB can be recovered, of course. After this result, 
we believe that obtaining gauge theories for non-linear systems is now an easier task. 

It is obvious that other and different definitions could be used with the same objective. 
For example, the generalized Euler formula, Abel or Fourier integral representation, Sonin, 
Letnikov, Laurent, Nekrasove and Nishimoto representation can be used. 

We presented an example where a version of the free relativistic particle was con- 
sidered. We calculated the Dirac bracket in the fractional embedding context and its form 
is very reasonable considering the conditions imposed, it is easy to see that by the first 
prescription we can not pinpoint all the physical features of the results. Having said that, 
we deem that the fractional embedding is a more general prescription encompassing even 
the approach introduced here. 

We strongly believe that quantization in a fractional context is an open area and deserves 
quite attention. We do not know yet the whole kind of problems that can be handled using 
this approach. Research in gravitation, condensed matter and field theory seem to be ready 
to be reinterpreted using the formalism of fractional calculus. 
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So far, it is not well known how to deal with dissipative systems. There are many paths of investi- 
gation in the literature and none of them present a systematic and general procedure to tackle the 
problem. On the other hand, it is well known that the fractional formalism is a powerful alternative 
when treating dissipative problems. In this paper we propose a detailed way of attacking the issue 
using fractional calculus to construct an extension of the Dirac brackets in order to carry out the 
quantization of nonconservative theories through the standard canonical way. We believe that us- 
ing the extended Dirac bracket definition it will be possible to analyze more deeply gauge theories 
starting with second-class systems. 
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I. INTRODUCTION 

The Dirac approach was very popular in the nineties, 
when an industrial production of papers concerning 
methods treating constrained systems were developed. 
The Dirac brackets (DB) [l| were an unmodified com- 
mon point between all papers about the subject. The 
motivation of many works were to convert second-class 
systems into first-class ones. The main objective was to 
obtain a gauge theory (first-class system), the holy grail 
for the Standard Model. Although not so popular as 
before, the analysis of constrained systems still deserves 
some recent attentions in the literature 0, Q . 

In few words we can say that the main feature of gauge 
theories is the existence of constraints which fix bound- 
aries in the phase space of gauge invariant systems to a 
submanifold [3]. Dirac covered all the main issues con- 
cerning constraint systems [l| , namely, a Hamiltonian ap- 
proach to gauge theories and general constrained theories 
and, consequently, the corresponding operator quantiza- 
tion procedure. Later on, the path integral method was 
found to be useful for quantizing gauge theories 043 and 
so-called second-class systems [8, 9], where the conven- 
tional Poisson bracket must be replaced by the DB in the 
quantization procedure. 

However, in constrained systems, it is possible to solve 
constraint equations The formalism proposed by 

Dirac for classical second-class constrained systems uses 
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the DB to deal with the evolution problem. The proce- 
dure is to apply the DB to functions of canonical variables 
in the unconstrained phase space, which avoid problems 
concerning the restriction of systems to constraints sub- 
manifolds 

On the other hand, there are various problems when 
considering classical systems besides the ones involving 
the quantization of second-class systems as we just have 
seen above. These problems encompass nonconservative 
systems. The curiosity about them is that the great ma- 
jority of actual classical systems is nonconservative and 
nevertheless, the most advanced formalisms of classical 
mechanics deals only with conservative systems [lo| . 

Dissipation for example, is present even at the micro- 
scopic level. There is dissipation in every non-equilibrium 
or fluctuating process, including dissipative tunneling 
[TH [T^ , electromagnetic cavity radiation [ij, [T^] and so 
on. 

One way to treat adequately nonconservative systems 
is through fractional calculus (FC) since it can be shown 
that, for example, a friction force has its form resulting 
from a Lagrangian containing a term proportional to the 
fractional derivative, which is a derivative of any non- 
integer order [loj . 

Non-linear dynamics is today an important subject of 
study in different physical and mathematical disciplines. 
However, its real success and a radically new understand- 
ing of non-linear processes occurred in the last 40 years. 
This understanding was inspired by the discovery and in- 
sight of a new phenomenon known as dynamical chaos. 
The reason for that is easy to understand, since any typ- 
ical system with more than one degree of freedom pos- 
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sesses chaotic motion for some initial conditions. We still 
do not know what is the measure of chaotic trajectories, 
but it seems that it is non-zero, and that makes the study 
of chaos important for constructing models of dynamical 
processes in nature [15]. 

FC is one of the generalizations of classical calculus. 
It has been used in several fields of science. FC provides 
a redefinition of the mathematical tools and it seems 
very useful to deal with anomalous and frictional sys- 
tems. In particular we can cite the continuous time 
random walk scheme as a physical counterpart exam- 
ple, where within the fractional approach it is possible 
to include external "fields" in a straightforward man- 
ner. Also the consideration of transport in the phase 
space spanned by both position and velocity coordinates 
is possible within the same approach. Moreover, the cal- 
culation of boundary value problems is analogous to the 
procedure for the corresponding standard equations fl^- 
|21|. Other important applications can be found investi- 
gating response functions where many studies have been 
reported on the phenomenon of nonexponential, power- 
law relaxation which is typically observed in complex sys- 
tems such as dielectrics, ferroelectrics, polymers and so 
on. The main feature of such systems is a strong (in 
general, random) interaction between their components 
in the passage to a state of equilibrium. The FC ap- 
proach to describing dynamic processes in disordered or 
complex systems such as relaxation or dielectric behav- 
ior in polymers or photo-bleaching recovery in biologic 
membranes has proved to be an extraordinarily successful 
tool. Some authors have proposed some fractional relax- 
ation models to filled polymer networks and investigate 
the dependence of the decisive occurring parameters on 
the filler content [H, [2^ . The study of exactly solvable 
fractional models of linear viscoelastic behavior is an- 
other successful field of application. In recent years both 
phenomenological and molecular-based theories for the 
study of polymers and other viscoelastic materials came 
up with integral or differential equations of fractional or- 
der. Some current models of viscoelasticity based on FC 
are usually derived from the Maxwell model replacing 
the first order derivative d/dt by its fractional version 
d" /dt", where a is not integer 124]. 

In this work we will use the well-known FC to analyze 
the well-established DB. The objective is to construct a 
generalized DB capable of treating a bigger number of 
mechanical systems than the standard DB. 

Since we believe that the FC has not been explored 
enough in field theory research yet, we tried to construct 
a self-sustained paper so that the issues are distributed as 
follows. In Section II we furnish a short history about FC 
together with its main equations and formulations. We 
will follow here the Riemann-Liouville (RL) approach. 
In Section III we establish the so-called fractional vari- 
ational principle, the ground stone for our cherished re- 
sult. However, we have to perform a modification of this 
fractional principle in order to include constrained sys- 
tems. In Section IV we analyze the same question but 



considering the action functional with generalized coordi- 
nates embedded in a fractional context, then we reinter- 
pret these different initial conditions to obtain a general 
formulation for the Dirac description for constrained sys- 
tems. In Section V, we use the free relativistic model to 
apply the fractional bracket. As usual, the conclusions, 
perspectives and last comments are depicted in the last 
Section, the sixth one. 



II. MODIFIED VARIATIONAL PRINCIPLE 

Nowadays, the inter- and multidisciplinarity among 
areas must be ever considered, and it therefore can be 
quite useful to study several problems from different ar- 
eas of science (besides the ones mentioned in the last 
section) such as viscoelasticity and damping, glassy con- 
densation, diffusion and wave propagation, electromag- 
netism, chaos and fractals, heat transfer, biology, elec- 
tronics, signal processing, robotics, system identification, 
genetic algorithms, percolation, modeling and identifica- 
tion, telecommunications, chemistry, irreversibility, con- 
trol systems as well as engineering, economics and finance 

Mm- 

It is well-known too from the current literature that the 
fractional approach can describe more precisely a myriad 
of physical systems. The formalism can be incorporated 
in many classical and quantum systems as described in 
the last section. We believe that its use can be extended 
up to field theory domain. 

The generalization of the concept of derivative with 
non-integer values goes back to the beginning of the the- 
ory of differential calculus. Nevertheless, the develop- 
ment of the theory of FC is due to contributions of many 
mathematicians such as Euler, Liouville, Riemann, and 
Letnikov ^t^. 

Since 1931, when Bauer 311 showed that we cannot 



use the variational principle to obtain a single linear dis- 
sipative equation of motion with constant coefficients, a 
new horizon of possibilities has been glimpsed. Nowa- 
days it has been observed that in physics and mathemat- 
ics the methodology necessary to understand new ques- 
tions has changed towards more compact notations and 
powerful nonlinear and qualitative methods. Derivatives 
and integrals of fractional order have been used to un- 
derstand many physical applications. For instance, ques- 
tions about viscoelasticity and diffusion process may have 
a more detailed description when this approach is used. 
In nature, the majority of systems contains an internal 
damping process and the traditional approach based on 
energy aspects cannot be used everywhere to obtain the 
right equations of motion. 

So, after Bauer's corollary, Bateman [32] proposed 
a procedure where multiple equations were obtained 
through a Lagrangian. Riewe [Ifll observed that using 
FC it was possible to obtain a formalism which could be 
used to describe both conservative and nonconservative 
systems. Namely, using this approach one can obtain the 
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Lagrangian and Hamiltonian equations of motion also for 
nonconservative systems. Agrawal studied a fractional 
variational problem [sl]- A fractal co ncep t applied to 
quantum physics has been investigated |34l |. 

The solution of a fractional Dirac equation of order 
a = 2/3 has been introduced [s^ and recently this sub- 
ject has been revisited [36| 



A. The fractional calculus 

We give a short introduction to FC. We believe that 
it will not provide the interested reader with all the FC 
tools, but we want to explain at least what is a fractional 
derivative. 

The first way to formally introduce fractional deriva- 
tives proceeds from the repeated differentiation of an in- 
tegral power 



(1) 



rfx" (to — n) ! 

For an arbitrary power fi, repeated differentiation gives 



r(A^ + i) 

T{fi-n+l) 



(2) 



with gamma functions replacing the factorial. The 
gamma functions allow for a generalization to an arbi- 
trary order of differentiation a, 



r(M + i) 

T{fi-a + l) 



(3) 



Of course, the objective of this work is that a can be a 
real number. The extension defined by the latter equa- 
tion corresponds to the RL derivative. It is sufficient for 
handling functions that can be expanded in Taylor series. 
A second way to introduce fractional derivatives uses the 
fact that the n-th derivative is the inverse operation to an 
n-fold repeated integration. Basic is the integral identity 

dy„ ■ ..dyifiyn) 

\yf{y){x~yr-K (4) 



{n-l)\ A 

A generalization of the expression allows one to define a 
fractional integral of arbitrary order a via 

aD-'^.f{x) = -i^ f dyf{y){x-yr'^ , (x > a). 

(5) 

A fractional derivative of an arbitrary order is defined 
through fractional integration and successive ordinary 
differentiation. 

For the time being we will not use the a-D" notation 
to indicate a fractional derivative. We will return to this 
in the future. 



After these few words about the fractional formalism, 
we think that it is important to justify our choice of using 
the RL fractional derivative instead of other very popu- 
lar fractional derivatives, as the Caputo time derivative, 
for instance. This last one could be another way to pur- 
suit the construction of the fractional DB. However, we 
understood that Caputo is more appropriated to appli- 
cations in several engineering problems due to the fact 
that it has a better relation with Laplace transform. 

In the last years the Caputo approach has been favored 
relatively to the RL one, because it is believed that the 
RL case leads to initial conditions without physical mean- 
ing. This was contradicted by Heymans and Podlubny 
37| that studied several cases and gave physical meaning 
to the RL initial conditions 'ss']. 

Still trying to clarify our objectives, let us affirm that 
Caputo's definition considers the differentiation inside 
the integral in order to solve the constant derivative prob- 
lem. For us it is an inconvenient way to describe some 
Lagrangian systems and consequently to obtain a con- 
sistent definition for our Dirac bracket, since we will be 
dealing only with field derivatives. For our main purpose, 
the RL approach is more convenient than Caputo one. 

However, it is well-known that several definitions of 
fractional derivatives and integrals exist. For instance, 
Griinwald-Letnikov, Caputo, Weyl, Feller, Erdelyi-Kober 
and Ricsz fractional derivatives as well as fractional Li- 
ouvillc operators, which have been popularized when 
fractional integration is performed in dynamical systems 
^,i4Q]. There is no equally simple definition that applies 
both to functions expressed as exponentials and to func- 
tions expressed as powers. In order to obtain a definition 
that is as general as possible, in order to be possible to 
attack other problems too, it has become conventional 
to use an integral representation discovered by Li- 
ouville 41 1 and extended by Riemann [l^. This is the 
main reason of our choice among other formulations of 
FC. Since we consider this work as the first step in the 
direction of the analysis of quantum field theories, we 
believe that a general definition of the fractional deriva- 
tive is the more convenient one. We will talk more about 
these perspectives in the last section. 



B. Modified Euler-Lagrange Equations 

Let us consider a smooth Lagrangian function. For 
any smooth path q : [a,b] — >■ M satisfying boundary con- 
ditions q{a) = Qa and q{b) = qh, consider an action- like 
Riemann Liouville fractional integral as considered in Eq. 
dSD above [H, 



1 



r(«) 



L{q{T),q{T),T){t- 



(6) 



where T(a) is the Euler gamma function, with a € (0, 1] 
and q — ^ is the derivative with respect to the intrin- 
sic time r G {a,t') and t € [to,t'] is the time for some 
observer in a particular reference system. 
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Notice that the Lagrangian in Eq. ([6]) is weighted 
by r^(^ — t)"~^. The time weighting acts as a time- 
dependent damping factor [s^ , and obviously when a — >■ 
1 we re-obtain the usual functional 



aSr[q{T)] = / L{q{T),q{T),T)dT. 



(7) 



Constructing the variation of the action functional, 
SSa = 0, we obtain after standard calculus the Euler- 
Lagrange equations associated with the fractional action 
integral ^ in Eq. © that 



dL 



dr \dqi 



1-adL 
t-T dqi 



= 0, i = 1, • • • ,n . (8) 



The Euler-Lagrange equation above, for some fractional 
action functional, must be obeyed. Now we will consider 
the invariance in phase space 



, 



(9) 



because we are interested in Dirac's quantization ap- 
proach. We intend to explore the same idea when con- 
strained systems will be under consideration. 
We can write, from Eq. ([7]), that the variation is 



Sasr = 



so that 



r( 



^ 6 f [pq - H{p,q,T)]{t-Tr-'dT ^ 

(10) 



r(a) Ja 



, 



(11) 



where L — pq — H(p,q,T) and the rest of the calcula- 
tion is standard from the variational calculus textbooks. 
The modification is due to PC formalism. However, it is 
direct to deal with this additional factor. Hence, after 
performing the variation of the Lagrangian as in Eq. (3) 
and of the damping factor; and isolating the coefficients 
for Sq and 6p that will be equal to zero, we obtain a new 
set of perturbed equations of motion, 



dH 

dpi 
dH 
dqi 



■Pi, 



(12) 
(13) 



which can be understood as the (fractional) Hamilton- 
Jacobi equations when this new action functional is con- 
sidered. It is clear that when a — !■ 1 our results will turn 
back to the usual case. 

We will see later that the expression Pi ^'^^^ i™" 
portant in our fractional DB formulation. The order a 
will be directly related to the fractional approach. The 
presence of a fractional factor is responsible for the 
generation of a time-dependent damping into the dynam- 
ics of the system, which is very useful to study models 



with smooth turbulence. Furthermore it is possible to es- 
tablish a relationship between the fractional Rayleigh dis- 
sipation function and the Euler-Lagrange equation [sll , 

_ ^ /aL\ aR ^ ^ ^ (14) 

dqi dr \ dqi J dqi 

where R is the fractional Rayleigh dissipation function 
given by 



R 



1 - a 

t-T 



L . 



(15) 



Note that in Eq. the dissipation function is part 

of the extended Euler-Lagrange equation. The origin of 
the third term is non-standard and is due to fractional 
analysis. 



III. MODIFIED VARIATIONAL PRINCIPLE ON 
CONSTRAINED SYSTEMS 

Now our main objective is to obtain an extended analy- 
sis which allows the quantization of classical systems with 
turbulence flow in field theory. We know that the quan- 
tization of a classical field theory in a natural context is 
not a straightforward unique process. The replacement 
of classical Poisson brackets by commutators of quantum 
operators cannot be carried out simultaneously for all 
conceivable dynamical variables without paying a price, 
i.e., internal obstructions will occur ,43, 44]. 

In general, the commutation formalism is restricted 
firstly to a certain class of variables, such as the canoni- 
cal coordinates of the theory. All commutators obtained 
will be derived from this first set. However the classi- 
cal theory may be substantiated in terms of any set of 
canonically conjugated variables in such a manner that 
the transition from Poisson brackets to quantum commu- 
tators leads to a "weird" quantum theory, depending on 
the chosen canonical coordinates system. 

This kind of problem usually occurs when the classical 
theory has constraints, and the right prescription for this 
was first formulated by Dirac 43[ and Bergmann and 
Goldberg [i^l, where they pinpointed the right bracket 
algebra to be used. Thus our goal now is to extend our 
last result to constrained systems. The action in Eq. pH)) 
can be considered in phase space. 



aSf 



1 



Via) 



pq - Hip,q,(t>a) it-Tr~'dT (16) 



where H is 



H = H + Xa(j)a ■ 



(17) 



The question involved in such systems is that when we 
carried out the Legendre transformation where L(q,q,t) 
becomes H{p,q,t), and defined the canonical momenta 
as Pi = ^perhaps the TV quantities are not all inde- 
pendent functions of the velocities. We cannot eliminate 
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the gi's and obtain M constraints equations 0a (9, p) = 0. 
Extending our discussion, we write the variation for Eq. 
(fTell as 



Sasr 



r(a) 
1 

rR J a 

dH 



pq - H{p,q,(l)a) 



Spq- [p-p 



1-a 



op oq op 



t - T 

Sp 



'dT 



Sci ] - ^Sq 



it- 



\Q-1 



= 



dr 



(18) 



After some algebraic manipulations, some terms can be 
I 



isolated allowing us to write the Hamilton- Jacobi equa- 
tions for the fractional constrained case. 



op Op 



P 



dH 
dq 



A, 



Oq 



1 -a 

t~T 



Pi , 



(19a) 
(19b) 



and again we have in Eq. (I19b[) a second term represent- 
ing the fractional contribution. 



The Fractional Dirac Bracket 



Consider a dynamical variable Q[pi,qi,t]; and using 
Eqs. we obtain 



dB de 

dt dq. 



dp, 



■Pi 



56 
'dt 



dQ (dH^^ ^ 

dqi \ dp^ ° dp^ J dpi 
de dH do dH 
dqi dp^ dpi dq^ 



dq^ ° dq^ ^ t — T 



+ Aa 



96, 



de 



^{e,H} + Xa{&,cj,a}-Pi 



ji dp"- dpi 
a-ldQ de 



Pt 



96 

a-lde 



t — T dp^ 



t ~ T dp^ 



dt 



dQ 
'dt 



(20) 



The constraints are dynamical variables too. Then, 
substituting some of the constraints in Eq. (|20p . we have 
that 

fA u\ , \ IJ, A \ a-ld(j)a d(j)a 

— = {0a, H} + A,{0a>.} ~P^J^^ + ^^ 

(21) 

and solving for Aa , we finally obtain a new result for the 
DB in a fractional context. 



{F,G}* = {F,G} 



PB 



{F,MC-,'p 



_i a — 1 d(j)b a — 1 dF 



t-T dpi 



t-T dpi 



(22) 



Our calculations show precisely this new result as a natu- 
ral extension for the DB. We must observe that, the usual 
DB appears inside the fractional correction and the ma- 
trix Cab = {0a, 0f)} is the Constraint matrix. It is obvious 
that when a — )• 1 we re-obtain the usual approach. 



IV. FRACTIONAL EMBEDDING 

Our next step is to build a general way to obtain the 
Dirac description for constrained systems. For this we 
will consider the problem under different initial condi- 
tions. A different and more general approach to analyze 



any dynamical system begins by considering the action 
as a function of generalized coordinates (loj . 

5[g(r),g(T)] - f L{ql{T),Ql,{T),T)dT (23) 

J a 

q:,^{aD'^rxr{t), qi,[t) = wr'xrit), 

with r — 1,2,...R coordinates considered, n = 
1,2, .. .M is the sequential order of the derivatives for 
the generalized coordinates q and n' = 1,2,... M' is the 
same for the coordinates Q. It can be shown [4^ that 
the necessary condition for an extremum of S is satisfied 
by 



0, (24) 



n— 1 '^'"^ n' — l 

and the momenta have the form 

N 



Pn- 1^ itDb) g^^^, 



k=n+l 
N 

E 

k=7i' + l 



dQL- 



(25) 



It is important to observe that we could extend the 
approach to a phase space just considering the usual ac- 
tion functional depending on the generalized fractional 
coordinates. 
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The Dirac formalism can be easily obtained here. It 
is well-known that it is useful in Lagrangian constrained 
systems. Now we propose its extension using the FC 
to encompass constrained non-conservative systems. Of 
course we could define our initial conditions in a different 
way and consequently obtain other final expressions. We 
realize that it is a very general form to deal probably 
non-linear systems and other kinds of phenomena. With 
this objective we define our constrained Hamiltonian, 

7? = if + ^AfeO, + (26) 

k k' 

where now the constraints are in fractional form too. We 
can define them by means of the RL prescription, 

k rt 



1 



T{k -a) \dt 
1 



(<-r)'=-"-Vfc(p,g,T)dT(27a) 



X' = 

T{k'-a) V dt 



'^^ ^\t~Tf-^-^Xk.{^,Q.T)dT. 



(27b) 



The resulting action is 



dt 



R N-l 



R N'-l 



2^ 2^ VWn +2^ ^ 

r — 1 71—0 r — 1 n'— 



H 



(28) 

and using of the variational principle, 5S — Q again, we 
can calculate the Hamilton- Jacobi equations. 



tD^Qn — 
t,L>„ TT„, — 



dH 
dH 

dH 

dpi 



Vk' 



Vk 



dqn 
dXk- 

dXy 



dQl 
d^k 



(29) 



This form of the Hamilton- Jacobi equations is new in 
the literature and introduces an extension of the Pois- 
son bracket into the RL context presented in Ref. '45']. 
It is natural that the next step is to obtain the proper 
DB expression. One way to do that is to consider some 
dynamical variable F{q^,p^^,Q'[^, ,t:'^) where 



dt dq, 
dF 



dF 



dQl, 



DbQn + 



tF>a 7r„' 



dF 
'dt ' 



(30) 



and after using Eq. ([2^ it is straightforward to build 
our final and main result for the DB in the RL context, 
namely 



{A,Br = {A,B} - {A,4>k}C^i{<l>i,B} 
^{A,Xk'}E,:,],{xi',B}, 



(31) 



where C and E are constraint matrices as in the standard 
Dirac constraint formalism. The consequent quantization 
can be described also in the standard way as 



[A,B] = ih{A,B}* 



(32) 



Now that we have constructed a proper fractional form 
for the DB, we believe that the conversion methods for 
obtaining first-class systems from second-class ones with 
non-linear models can be carried out. 

In the next section we will apply our result obtained 
in Eq. ()3ip in a well-known and simple model, the rela- 
tivistic free particle. 



V. THE RELATIVISTIC FREE PARTICLE 

Our objective in this section is to study one applica- 
tion in the light of the fractional DB introduced in Eq. 
([3T|). To fix our ideas developed before, we will consider 
a simple example, the relativistic free particle model, to 
apply the fractional embedding. 

This model is well-known, and its usual Lagrangian is 
given by 



(33) 



Using the ideas of the last sections, the action under 
consideration is 



L 



-TO^[(ai?"t)"^r]' + [(tD"fa)"'x,]'. (34) 



We will restrict our calculations to the case when a = 
1/2 and n = n' = 1. Therefore we have 

L = -'I 
or in a simpler way 



1 1 r 1 

{aD2t)Xr + {tD2,,)xr 



(35) 



L = -m^{q\r + {Q\r ■ 



(36) 



Of course we could consider different orders for the 
derivative operator, but our main intention now is to ap- 
ply the method and to show its usefulnes to canonically 
quantize fractional systems. 

Using the definition of generalized momentum we ob- 
tain the two conjugated momenta. 



Po = -"^ 



Q\ 



The primary constraint is 



2,2 2 
Po -f TTq - TO 



(37a) 
(37b) 

(38) 
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As the canonical Hamiltonian is zero, to construct the ex- 
tended Hamiltonian by Dirac's prescription we can write 



and since it is a first-class constraint 



(39) 



(/-i = {01, i?} = 0; 

therefore we do need to fix the gauge and our choice is 

9? + 0!-r = 0. (40) 

Now we have two second-class constraints 

{01,02} - -2(pO + vrg), (41) 

and the new extended Hamiltonian can be written as 

H [ipir + {irlf - m^] + + {Q\f ~ r] . 

(42) 

The time evolution of these constraints give us the correct 
form of the Lagrange multipliers. The extended Hamil- 
tonian in its final form is 



H = 



1 



2(p1 



(43) 



Using the DB definition from Eq. ([3T|) we can calculate 
finally that 



{qi,Q\r = {K,^ir = 0, 

Po ^0 

Po - 



(44a) 
(44b) 

(44c) 



and the quantization is directly obtained using the stan- 
dard Eq. ([5^ . We can observe that the brackets ob- 
tained above have (a kind of) expected results. In other 
words, the commutative result in Eq. (I44al) is standard. 
The results in Eq. (j44bp and Eq. (|44cp are standard also 
in the first term. The second terms in both these equa- 
tions are consequences of the fractional approach. The 
result obtained in Eq. (|44a|) make us think about non- 
commutative issues. We will talk more about this in the 
next section. 

We have to clarify the interested reader that here we in- 
troduced the EC to investigate nonconservative physical 
systems. Consequently, we are beginning to fathom other 
physical features inside the fractional formalism, different 
from the current literature. To help us in this task, the 
next step would be to investigate a solid nonconservative 
system like radiation damping, which is not completely 
understood using the standard (non-fractional) formal- 
ism. This is a target of current research by the authors 
and will be published elsewhere. 



VI. CONCLUSIONS AND PERSPECTIVES 

Our main motivation is to develop an approach 
based on fractional variational calculus to handle non- 
conservative constrained systems, since EC can be used 
to deal with a frictional force, for example. In other 
words, our effort is to construct a fractional DB. Con- 
sequently, we showed that it is possible to think about 
quantization in this scenario. In this way we proposed 
two kinds of fractional formulations for the DB. 

The first one is based on the RL derivative, but it is 
incorporated directly into the action functional. We ob- 
tained the Hamilton- Jacobi equations that are deformed 
by the fractional contribution. Consequently the DB also 
has the same kind of modification. 

However, this first approach does not seem to be the 
right one. Therefore, we changed the formalism con- 
sidering a usual form for the action, but redefining the 
coordinates in a generalized prescription using the frac- 
tional definition according to Riewe's prescription. The 
constraints were defined in the same way and the con- 
sequence was the extension of the usual Euler-Lagrange 
equations of motion into a fractional scenario. 

We obtained the final form for the fractional DB which 
has an additional term due to the EC contribution. We 
showed that the standard DB can be recovered, of course. 
After this result, we believe that obtaining gauge theories 
for non-linear systems is now an easier task. 

Other and different definitions could be used with the 
same objective. Eor example, the generalized Euler for- 
mula, Abel or Fourier integral representation, Sonin, Let- 
nikov, Laurent, Nekrasov and Nishimoto representation 
can be used. 

To apply our fractional Dirac formalism, we used an 
example where a version of the relativistic free par- 
ticle was considered. We calculated the respective DB 
in the fractional embedding context, and its form results 
very reasonable considering the conditions imposed. It is 
reasonable to see that with this first prescription for the 
fractional DB we cannot pinpoint all the physical features 
of the results. 

Geometrically speaking, fractional models provide us 
with a memory effect about the convolution integrals and 
give us some differential equations with a bigger expres- 
sive power. This allows us to consider several different 
physical situations such as viscoelasticity and more ab- 
stract scenarios such as mapping using tensorial fields. 
Physically we can understand the derivative of order a 
of some individual velocity v as the same velocity Voh 
from the point of view of an independent observer 46] . 

As a perspective, one possible target of research would 
be the noncommutative fractional dynamics. We can ob- 
serve from Eq. (|44ap that this algebra is commutative. 
However, we can ask whether this is the standard pattern 
or whether there exists any algebra induced by EC where 
something like Eq. (|44ap is not true, like in string theory 
using ordinary calculus) in a magnetic field background 
47| . Considering this scenario, we can ask whether the 
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application of the symplectic formalism (where we can 
choose the zero-mode, like in Ref.psj) coupled with frac- 
tional formalism, the canonical noncommutativity can 
be obtained. This is our next task, together with the 
construction of a Moyal-Weyl product using FC to han- 
dle problems such as noncommutative quantum mechan- 
ics where the Moyal-Weyl product will have a fractional 
form and maybe we can obtain interesting results com- 
paring with the fractional quantum mechanics ones in 
the current literature. Or with a noncommutative alge- 
bra obtained with the DB obtained here. One manner 
that we have to re-obtain the commutativity introduces 
a noncommutative parameter in the original non-linear 
system. We can investigate how it works in a fractional 
scenario. 



Besides the applications proposed here, we strongly be- 
lieve that quantization in a fractional context is an open 
area and deserves more attention. We do not know yet 
the whole kind of problems that can be handled using 
this approach. In general, research in gravitation, con- 
densed matter and field theory seem to be ready to be 
reinterpreted using the formalism of FC. 



VII. ACKNOWLEDGMENT 

The authors would like to thank the referees. Their 
suggestions and concerns contributed substantially to 
transform our paper into a better one. 



[1] P. A. M. Dirac, "Lectures on Quantum Mechanics" (Bef- 
fer Graduate School of Science, Yeshiva University, New 
York, 1964). 

[2] T. G. Philbin, New J. Phys. 12, 123008 (2010). 

[3] M. Maeno, Prog. Theor. Phys. Suppl. 188, 217 (2011). 

[4] M. I. Krivoruchenko, A. A. Raduta and A. Fraessler, 
Phys. Rev. D 73, 025008 (2006). 

[5] M. Henneaux, Phys. Rept. 126, 1 (1985). 

[6] L. D. Faddeev and A. A. Slavnov, "Gauge fields, In- 
troduction to Quantum Tfieory" (Addison- Wesley, New 
York, 1988). 

[7] M. Henneaux and C. Teitelboim, "Quantization of Con- 
strained Systems'" (Princeton University Press, Prince- 
ton, New Jersey, 1992). 

[8] I. A. Batalin and E. S. Fradkin, Phys. Lett. B 180, 157 

(1986) . 

[9] I. A. Batalin and E. S. Fradkin, Nucl. Phys. B 279, 514 

(1987) . 

[10] F. Riewe, Phys. Rev. E 53, 1890 (1996); Phys. Rev. E 

55, 3581 (1997). 
[11] A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 

211 (1981); Ann. Phys. 149, 374 (1983). 
[12] J. Ankerhold, H. Grabert and G. L. Ingold, Phys. Rev. 

E 51, 4267 (1995). 
[13] L R. Senitzky, Phys. Rev. 119, 670 (1960). 
[14] K. W. H. Stevens, Proc. Phys. Soc. London 72, 1027 

(1958). 

[15] G. Zaslavsky, "Hamiltoman Chaos and Fractional Dy- 
namics" (Oxford University Press, New York, 2005). 

[16] R. Metzler, E. Barkai and J. Klafter, Phys. Rev. Lett. 
82, 3563 (1999). 

[17] R. Metzler and J. Klafter, Phys. Rept. 339, 1 (2000). 

[18] R. Hilfer, Chem. Phys. 284, 399 (2002). 

[19] R. Metzler and J. Klafter, J. Phys. A 37, R161 (2004). 

[20] D. Fulger, E. Scalas and G. Germano, Phys. Rev. E 77, 
021122 (2008). 

[21] M. A. E. Herzallah and D. Baleanu, Nonlinear Dyn. 58, 
385 (2009). 

[22] R. Metzler, W. Schick, H. G. Kilian and T. F. Nonnen- 
macher, J. Chem. Phys. 103, 7180 (1995). 

[23] A. Stanislavsky 1, K. Weron and J. Trzmiel, Europhys. 
Lett. 91, 40003 (2010). 



[24] W. G. Glockle and T. F. Nonnenmacher, Macromolecules 
24, 6426 (1991). 

[25] J. A. Tenreiro Machado, M. F. Silva, R. S. Barbosa, 
I. S. Jesus, C. M. Reis, M. G. Marcos and A. F. Gal- 
hano, "Some Applications of Fractional Calculus m Engi- 
neering" , Mathematical Problems in Engineering, 639801 

(2010) . 

[26] J. A. Tenreiro Machado, V. Kiryakova and F. Mainardi, 
Commun. Nonlinear Sci. Numer. Simulat. 16, 1140 

(2011) . 

[27] K. B. Oldham and J. Spanier, "The Fractional Calculus: 
Theory and Applications, Differentiation and Integration 
to Arbitrary Order" (Academic Press, New York, 1974). 

[28] K. S. Miller and B. Ross, "An Introduction to the Frac- 
tional Calculus and Fractional Differential Equations" 
(Wiley-Interscience, New York, 1993). 

[29] I. Podlubny, "Fractional Differential Equations. An In- 
troduction to Fractional Derivatives, Fractional Differ- 
ential Equations, Some Methods of Their Solution and 
Some of Their Applications" (Academic Press, San 
Diego, 1999). 

[30] A. Kilbas, H. Srivastava and J. Trujillo, "Theory and Ap- 
plications of Fractional Differential Equations" (Elsevier, 
Amsterdam, 2006). 

[31] P. S. Bauer, Proc. Natl. Acad. Sci. USA 17, 311 (1931). 

[32] H. Bateman, Phys. Rev. 38, 815 (1931). 

[33] O. P. Agrawal, J. Math. Annal. Appl. 272, 368 (2002). 

[34] N. Laskin, Phys. Lett. A 268(3), (2000) 298; Chaos 10, 
780 (2000). 

[35] A. Raspini, Phys. Scripta 64, 20 (2001). 
[36] D. W. Dreisigmeyer, P. M. Young, J. Phys. A 36, 829 
(2003). 

[37] N. Heymans and L Podlubny, Rheol. Acta 37, 1 (2005). 

[38] M. Ortigueira and F. Coito, "The Initial Conditions of 
Riemann-Liouville and Caputo Derivatives" (in the 6th 
EUROMECH Conference ENOC, Saint Petersburg, Rus- 
sia, 2008). 

[39] R. A. El-Nabulsi, D. F. M. Torres, J. Math. Phys. 49, 
053521 (2008). 

[40] R. A. El-Nabulsi, Chaos Solitons Fract. 42, 2614 (2009). 
[41] J. Liouville, J. Ec. Polytech. 13, 1 (1832); 13, 71 (1832); 
13, 163 (1832). 



Fractional Dirac Bracket and Quantization for Constrained Systems 



9 



[42] B. Riemann, "Gesammelte Werke" (Teubner, Leipzig, p. 
353, 1876). 

[43] P. A. M. Dirac, Can. J. Math. 2, 129 (1950). 
[44] P. G. Bergmann and I. Goldberg, Phys. Rev. 98, 531 
(1955). 

[45] D. Baleanu and S. I. Muslih, in "Advances in Fractional 
Calculus, Theoretical Developments and Applications in 



Physics and Enginnering," edited by J. Sabatier et al. 

(Springer, Netherlands, 2007), p. 115. 
[46] I. Podlubny Fract. Gale. Appl. Anal. 5, 367 (2002). 
[47] N. Seiberg and E. Witten, J. High Energy Phys. 9, 32 

(1999). 

[48] N. Braga and C. F. L. Godinho, Phys. Rev. D 65,085030 

(2002). 



